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Abstract. A technique for determining the permeability of a phospholipid 
membrane on a single giant unilamellar vesicle (GUV) is described, which 
complements the existing methods utilizing cither a planar black lipid membrane 
or sub-micrometre-sized liposomes. A single GUV is transferred using a 
micropipette from a solution of a nonpermeable solute into an iso-osmolar solution 
of a solute with a higher membrane permeability. Osmotical swelling of the vesicle 
is monitored with a CCD camera mounted on a phase contrast microscope, and a 
sequence of images is obtained. On each image, the points on the vesicle contour 
are determined using Sobel filtering with adaptive binarization threshold, and 
from these, the vesicle radius is calculated with a great accuracy. From the 
time-dependence of the vesicle radius, the membrane permeability is obtained. 
Using a test set of data, the method provided a consistent estimate of the POPC 
membrane permeability for glycerol, P = 1.7 X 10~ 8 m/s, with individual samples 
ranging from 1.61 X 10 — 8 m/s to 1.98 X 10 — 8 m/s. This value is ta 40% lower 
than the one obtained on similar systems. Possible causes for this discrepancy are 
discussed. 
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1. Introduction 

All living cells are enclosed by a lipid bilayer, which serves as a barrier for 
macromolecules and ions. Small uncharged molecules can, however, permeate across 
the membrane without requiring assistance of transmembrane proteins. Permeability 
of a lipid bilayer has traditionally been measured either electrically on a planar 
black lipid membrane pQ, or optically on sub- micrometre-sized liposomes (e.g., large 
unilamellar vesicles, LUVs), through osmotic swelling or shrinking of vesicles and a 
subsequent change in the light scattering [2j [3] . While these techniques have proven 
to be successful, they have two drawbacks. Firstly, they offer only indirect insight into 
the mechanism involved, and secondly, an open question remains whether the results 
obtained on the membranes of sub-micrometre-sized LUVs can be applied onto the 
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membranes of micrometre-sized giant unilamellar vesicles (GUVs), i.e., whether the 
membrane permeability is influenced by its curvature. 

In this paper, we propose a complementary technique for measuring the 
permeability of a lipid bilayer on GUVs using an analysis of a sequence of 
videomicroscopy images. Processing and analysis of digital video images acquired 
by an optical microscope has become an indispensable tool in biophysics, in particular 
for the analysis of the shape of contours of phospholipid vesicles, either free [H [5] or 
aspirated in a micropipette [6] . 

2. Description 

2.1. Experimental setup for image acquisition 

In the experiment, a single spherical GUV made from l-palmitoyl-2-oleoyl-sn-glycero- 
3-phosphocholine (POPC) was selected, fully aspirated into a glass micropipette with 
a diameter exceeding the diameter of the vesicle, and transferred from a solution of 
a solute with a very low membrane permeability {e.g., glucose or sucrose) into an 
iso-osmolar solution of a solute with a higher membrane permeability {e.g., glycerol), 
where the content of the micropipette was released, and the micropipette was removed. 
Vesicle behaviour is recorded with a CCD camera mounted on the microscope. Upon 
the transfer, the vesicle swells osmotically until it bursts, whereupon another cycle 
of swelling is commenced. The rate of swelling depends on membrane permeability 
and the difference of the concentration of the permeable solution outside and inside 
the vesicle. Since the whole sequence of the vesicle response is recorded, and the 
vesicle radius can be determined from micrograph images with a great precision, an 
estimate for the concentration difference can be obtained from the vesicle geometry, 
thus making an estimate of the membrane permeability possible. Vesicle transfer 
from a solution of a non-permeable into a solution of a permeable solute rather than 
vice versa, thus resulting in osmotic swelling rather than osmotic shrinking, has been 
chosen since it completely determines the geometry of a vesicle, i.e., the vesicle is 
spherical throughout the experiment, thus making a high-precision measurement of 
its radius possible. The micromanipulation experiment is explained in detail in [7]. 

The experimental setup consisted of a phase contrast inverted optical 
microscope (Nikon Diaphot 200, objective 20/0.40 Ph2 DL; Tokyo, Japan) with the 
micromanipulating equipment (Narishige MMN-l/MMO-202; Tokyo, Japan) and a 
cooled CCD camera (Hamamatsu ORCA-ER; C4742-95-12ERG; Hamamatsu, Japan), 
connected via IEEE-1394 to a PC running Hamamatsu Wasabi software. The camera 
provides 1344 x 1024 12-bit grayscale images. In the streaming mode, the camera 
throughput is 8.9 images/s. Even though custom solutions have been shown to provide 
a higher throughput in similar cases [8], the software supplied by the manufacturer 
was proved to be adequate for the camera. 

2. 2. The model for vesicle inflation and burst 

The basic ideas of the model explaining the swelling-burst cycle will be presented first. 
After the transfer, a spherical vesicle filled with an impermeable solute is immersed in a 
solution of a permeable solute. Since the concentration of the permeable solute outside 
the vesicle (c p o) exceeds its concentration inside the vesicle (c p ), it tends to diffuse into 
the vesicle. During this phase, the amount of the impermeable solute inside the vesicle 
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Figure 1. A schematic indication of the vesicle radius R as a function of time 
t upon transfer from a solution of an impermeable solute into an iso-osmolar 
solution of a permeable solute. 



(N[) remains constant, while the amount of the permeable solute (N p ) increases. Its 
influx is accompanied by an influx of water required to maintain the osmotic balance. 
Based on the known data pQ , it is assumed in this model that the permeability of the 
membrane for water exceeds all other permeabilities by several orders of magnitude, 
and consequently the osmotic balance is achieved almost instantly. Consequently, in 
order to maintain the osmotic balance, the vesicle volume increases as well. In an 
instantaneous event of a vesicle burst, a fraction of the solution is ejected from the 
vesicle interior. During the burst, the concentration of both solutes (c p , c{) inside the 
vesicle stays constant, while their amount decreases, the decrease being proportional 
to the decrease of the vesicle volume. The cycle then repeats itself with new values 
of the concentrations c p and q. The theory of the osmotic inflation-burst cycle is 
described in more detail in [7] . 

The flux of the permeable solute into the vesicle (j) is proportional to the 
difference of the concentration of permeable solute outside (c pu ) and inside (c p ) the 
vesicle, 

j = P(c p0 - c p ) , (1) 

where P is the permeability of the membrane for the permeable solute. The number 
of its molecules inside the vesicle (N p ) increases with time at a rate: 



dN p 
dt 



PA 



■pO 



iVp 
V 



(2) 



Both the vesicle volume V and the membrane area A vary with time, which means 
that in general, |J2J) can only be solved numerically. 

For realistic parameter values, however, numerical integration is not needed. It is 
known [11] that the lipid membrane can only expand by approximately 4% before 
its tensile strength is reached ("stretched" radius R s ), at which point the vesicle 
bursts and ejects its excess volume, and its radius R returns to its "relaxed" value R r 
(figure [TJ. 

Such small changes of R can be linearized, leading to a linear dependence of V in 
the time between bursts. A rearrangement of Eq. 11 in [7] yields a relation between 
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the rate of the volume increase in the n-th cycle, AV n / At n , and the concentration of 

(n) 

the impermeable solute in the vesicle during this cycle, iVj : 
AK V s . n - V t . n 3P iV ; (n) 



At n — in -Rr,0 c p0 

We have denoted V^ )7l = 47ri2j? n /3 and I4„ = 47ri?^„/3. 

While t n , i? rj „ and i? SjT1 can be determined directly, the value of N> at the ra-th 
burst can be calculated. Since the amount of the impermeable solute inside the vesicle 
remains constant during the time between successive bursts, one can write: 

Af> = q <°V ri0 = c f»y s ,„ , 

= q (1 V tll = Cl (2 V M , (4) 
From here, a general expression for Cj can be obtained: 



n-1 

. (o) n ^ (5) 



c («) = c (o) TT }^ 



3=0 



The initial amount of the impermeable solute inside the vesicle is known, = 
V Tl oCt ■ We also know that the solutions inside and outside the vesicle are iso-osmolar 
(c| = c p o). Taking into account that the relaxed radii of the vesicle are approximately 
equal (V^ n rj K-,o), and substituting this into ([3]), one obtains: 

f ^ 4 <on^. (6) 

Thus, plotting Ay„/Ai„ against 4irRf YYjZo V Tj j/V at j, one should obtain a straight 
line with the slope equal to P. 

2. 3. Determining vesicle radius from a videomicrograph 

Phase contrast micrographs exhibit a distinct "halo" , if the index of refraction of the 
medium inside does not match the refractive index of the medium outside |12j . 

A Sobel edge detection operator [13] was used to determine the radius of the 
vesicle. The choice of threshold of the binarization applied to the image after the 
convolution affects the result in an important way. Since the contrast of the halo 
decreases throughout the time of recording, the binarization threshold needs to be 
adjusted dynamically in order to yield the desired result (figure [J). The binarization 
threshold was selected on the basis of the number of points it yields, in a way that 
the number of points making up the contour was between 2irR/3 and 2ttR 7 where R 
is the running average of previously determined values, expressed in pixels. In this 
interval, the calculated radius is virtually independent of the number of points taken 
into calculation (figure [3]). 

A circle with a centre (x c , y c ) and a radius R is fitted to a set of n points (xi, y{) 
obtained in the previous step. A computationally efficient algorithm for fitting a 
general conic to a set of data points [2] was adapted specifically for fitting ellipses 
[T5] . It turns out that the system obtained in [TS] is unstable, e.g., the approach does 
not work for the data, sampled from a circle. A more stable and efficient modification 
has been proposed by Hah'f and Flusser 16J. 
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(a) (b) 

Figure 2. Adaptive Sobcl threshold. As the composition of the solution in the 
vesicle interior approaches the composition of the vesicle exterior, the contrast 
fades from vidcomicrograph frames taken at the beginning of the experiment (a, 
top row) to frames taken towards the end of the experiment (b, top row); the 
second row shows the intensity profile of the vesicle cross-section shown in the 
same scale. It is thus necessary to dynamically adjust the binarization threshold 
of the Sobel edge operator. High values of threshold, which produce an acceptable 
image at the beginning of the experiment (a, third row), fail to discriminate the 
edge towards the end of the experiment (b, third row). On the other hand, low 
values of threshold, which are appropriate for the image frames towards the end 
of the experiment (b, bottom row) produce stray points at the beginning of the 
experiment. The algorithm for the dynamical adjustment of the binarization 
threshold selects an image corresponding to the one in the third row for (a), and 
an image corresponding to the one in the fourth row for (b). The bar in the top 
left frame represents 20 /im. 
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Figure 3. For a given image, the calculated radius of a vesicle depends on the 
number of points taken into consideration. When the binarization threshold is 
very high, very few points are selected and the calculated radius is significantly 
affected by any new point taken into consideration. In an intermediate region, 
the calculated radius does not depend on the number of points. Above a certain 
threshold value, spurious points taken into consideration affect the calculated 
radius. The data plotted are for the image plotted in figure [2^,; the interval 
[2ttR/3, 2ttR] approximately translates to [150,460], corresponding to the interval 
[9.90, 22.60] for the binarization threshold. 



Further simplification can be employed for fitting circles, which reduces the 
problem to solving a linear 3x3 system rather than solving a generalized eigenvalue 
problem, as required in [15) . If the matrix A and the vector b are constructed as 
follows: 



b = 



E* x iVi 

E* Xi 
E* 



E* aw 

E^ 2 
E*^ 

2 



Ei x i 

n 



vi 



(7) 



(8) 



then the required parameters are stored in a vector v = (v-i, 1)2, 1*3), obtained as the 
solution of a set of linear equations, 

A v = b . (9) 



In particular, 



Vi 



Vc 

R 



2 



'•■2 



(10) 

(11) 

(12) 



VZ + v 3 ■ 

The matrix A in |7]) is positive definite, provided that at least 3 data points 
are non-collinear. This can be easily seen from the factorization A = C T C, where 
C := [x y 1] with vectors of data points x, y, and a vector of ones 1. Note that the set 
of linear equations ©, even though derived by a different procedure, is equivalent to 
the one obtained by Kasa [17] by a modified least square criterion, i.e., by minimizing 
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Figure 4. The determined vesicle radius R as a function of time t upon a transfer 
from a 0.2 mol/L sucrose/glucose solution into an iso-osmolar glycerol solution. 

the sum J2ii( x i ~ x c) 2 + (y% — J/c) 2 — R 2 ] 2 ■ We would like to remark that a great deal 
of research has been done on the approximation of circular arcs and whole circles (see, 
e.g., 18] and the references therein). Since in our case the vesicle contour is known 
to have circular shape, special polynomial or spline approximations are not needed. 

A script was written in GNU Octave ( |http : //www . octave . or g/, [19]). a 
Matlab-compatible software package, to determine the radius of the vesicle from 
a videomicrograph image. Figure 2] shows a series of calculated radii from an 18- 
minute recording of a vesicle. As verified by a visual inspection, occasional apparent 
transient increases or decreases of the radius are due to image defocussing (as the 
cover slide, along which vesicles move, is not perfectly aligned with the focal plane, 
the experimentalist needs to occasionally adjust the focussing knob on the microscope 
in order to assure that the vesicle is still in focus). Such transient changes in the 
apparent calculated vesicle radius occur on an a time-scale of seconds and it is thus 
easy to discriminate them agains real changes of vesicle radius. 

Another GNU Octave script was used to read the output of the first script, i.e., 
the vesicle radius as a function of time, and determine the points at which vesicle bursts 
occurred. The algorithm employed was a simple one: the script looked for negative 
jumps in vesicle radius, the magnitude of which exceeded some pre-determined value 
and which did not follow the previous jump too soon. Both parameters — the minimum 
amplitude and the minimum time interval between the jumps — were hand-tailored to 
fit the experimental data. Using this algorithm, 90-95% of the jumps were detected; 
the missing ones were added by hand upon comparing the results to the data. 

3. Results and Discussion 

3.1. Test set of data 

The algorithm was tested on a set of data used in [7j (figure^ ■ Membrane permeability 
has been calculated for four recorded series of bursts (vesicles 2-5 in table 1, [7]), 
yielding an average value P = 1.7 x 10 -8 m/s at a room temperature (26 ± 2°C), with 
the estimates on individual vesicles ranging from 1.61 x 1(T 8 m/s to 1.98 x 1(T 8 m/s. 
This value is lower than the value (2.09 ± 0.82) x 10 -8 m/s obtained in [7J by another 




Figure 5. The rate of the volume increase AV n /Af n plotted against 
A-kR^ q n^Tt) 1 ^tj/^SiJ ' snown f° r four different vesicles transferred from 0.2 mol/L 
glucose/sucrose into 0.2 mol/L glycerol (denoted O (19 points), • (31 points), 
□ (17 points), and I, 20 points). Least-squares fitting of a straight line through 
the origin yields values of P (1.61 ± 0.07) X 10~ 8 m/s, (1.62 ± 0.04) X 10 -8 m/s, 
(1.98±0.07) X 10~ 8 m/s, and (1.93±0.08) X 10~ 8 m/s, respectively (plotted with 
a dotted line; the lower two lines are barely distinguishable). Here, the slope of a 
straight line through n points (xi,m) is calculated as m = XiVi / 53i , and 
its standard error as — mx i) 2 / {( n ~ 1) X D- A solid line fitted through 

all four sets of points corresponds to P = (1.69 ± 0.03) X 10 -8 m/s. 



method on the same data set, as well as from other published data, e.g., the value for 
DOPC, 2.75 x l(T 8 m/s at 30°C [3], and two values for egg-PC, (4.3 ±0.1) x l(T 8 m/s 
[2DJ and 5.4 x l(T 8 m/s at 25°C [21]. 

One has to note, however, that the only result against which we can directly 
compare the obtained result is the one obtained by a different analysis on the same 
data set [7]. Both results agree within the margins of standard error in [7 . All 
other results were obtained on systems which were different, even though they share 
some similarities. Paula et al. [5J did not include POPC in their measurement. 
The most similar system of those included is DOPC, which has two unsaturated 
acyl chains (POPC has one unsaturated and one saturated chain), ft is known 
that increasing chain saturation decreases membrane fluidity, and membrane fluidity 
has been shown to be well-correlated with the permeability of membrane for small 
uncharged solutes [22]. Furthermore, the experiment has been conducted at a 
slightly higher temperature, and it is known that membrane permeability scales with 
temperature approximately as Arrhenius factor, oc exp(—E a /kT) [23j . 

The other two measurements [20j [21] were conducted on natural egg 
phosphatidylcholine. According to the specifications provided by the producer (Avanti 
Polar Lipids, Inc.), this is a mixture of various lipids, both saturated and unsaturated, 
ranging in length of the acyl chain from 16 to 20 carbon atoms. It is to be expected that 
a heterogenous membrane has a higher permeability, and indeed, the permeabilities 
of egg-PC membrane listed in 20, 21] are consistently higher than those for DOPC 
membrane listed in [3]. 

As one can see, the presented method offers a possibility to determine membrane 
permeability with a relative accuracy comparable to or greater than other published 
results. Based on the above argumentation, one could expect that the value obtained 
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for permeability of POPC membrane should be slightly lower than the published 
values for other similar systems, and indeed it is. We have not shown, however, that 
the extent of the observed discrepancy can be explained solely by the fact that we were 
working with a different system. We can think of several possible additional factors 
which could also affect the estimated permeability. 

Firstly this discrepancy may indicate that the membrane permeability depends 
on the membrane curvature and thus one may indeed expect different values for the 
membrane of GUVs, which are essentially planar, than for the membranes of LUVs 
and other smaller aggregates. A comparison by Brunner et al. [53] however shows that 
planar membranes (black lipid membranes, BLM) exhibited higher permeabilities for 
ions and sugars than the membranes of sonicated small unilamellar vesicles (SUVs); 
it is to be noted, though, that the authors report problems with leakage and attribute 
a comparatively high value for BLM permeability to this artefact. 

Secondly, it may be due to the fact that the model presented in section 12.21 does 
not take into account the membrane elasticity. It has already been pointed out [55] 
that failing to account for membrane elasticity in the interpretation of the membrane 
permeability measurements in a similar experiment leads to an underestimate of 
membrane permeability. Our earlier estimates [7] however show that this correction 
amounts to only 2-3% and thus cannot explain the observed difference. 

Thirdly, the apparent value of permeability, as obtained in the experiment, is 
coupled with the diffusion of the permeable solute through the stagnant layer [26] 
and is therefore lower than the true permeability of the membrane alone. A quick 
estimate however shows that the effect of the stagnant layer only affects the calculated 
permeability to a few percent. 

Finally, one has to bear in mind that the phase contrast technique does not yield 
the true position of the vesicle membrane. Instead, we have to rely on the "halo" , 
which is an artefact of the technique, and define the position of the membrane as the 
point of the maximal gradient of the image intensity. For a straight-edge object, it has 
been calculated [27] that the maximal steepness of the profile slope indeed coincides 
with the position of the edge. A rigorous theoretical treatment would require repeating 
this calculation for each vesicle, taking into account the real parameters of the optical 
system used, which would considerably complicate the analysis without a promise of a 
significant improvement. As a consequence, even the most precise analyses of vesicle 
contours to date (28] [29] dismiss this discrepancy as negligible and rely to the result 
of Wilson and Sheppard [27] . All in all, we can conclude that the method presented 
in this paper needs to be applied to a wider variety of systems before this feature can 
be fully explained, and then possibly applied to other interesting phenomena, such as 
the critical phenomena predicted [30j . 

3.2. Comparison with other halo-based techniques 

Although it is based on previous works [101 [7] , the presented technique uses a different 
method for determining the composition of the vesicle interior, i.e., the amount of 
the impermeable solute JV; . The previously published methods relied on the height 
of the "halo" of the phase-contrast image for determining the composition of the 
vesicle interior and tacitly assumed that all the jumps of the vesicle radius were equal. 
Figure [JJ demonstrates that this is not necessarily the case. Furthermore, since no 
ab initio calculation of the "halo" profile has been done, the authors in [10] and [7] 
had to resort to the calibration of the halo. As the halo profile depends not only on 
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the refractive indexes of the solutions in the vesicle interior and the vesicle exterior, 
but also on the size of the vesicle and, in absolute terms, on the properties of the 
optical system used (e.g., & slight displacement of the phase ring from the optical 
axis affects it considerably), this analysis is cumbersome, and any laboratory trying 
to apply this technique needs to repeat the calibration on their system. By replacing 
the need for the measurements of the halo properties such as its height or width 
with a more precise measurements of the vesicle radius, which enables us to estimate 
the composition of the vesicle interior using simple geometry ((6]), we eliminated the 
need for the cumbersome calibration altogether, and thus made the technique much 
more portable. A further advantage is that shedding the reliance on the quantitative 
properties of the halo also opens the possibility for using techniques for improving 
uneven illumination (see, e.g., |31j). if necessary. Note, however, that when using the 
presented technique, it is important to record the vesicle behaviour right from the 
time of the transfer, as missing a burst yields to underestimating the permeability. 

Applying the technique from |10j to the measurement of small polar molecules, 
as has been done in [7J, introduces a further problem, since the calibration itself 
involves a tricky part — the experimentalist has to determine a suitable point for the 
initial concentration difference, i.e., & point late enough that the amount of the initial 
solution transferred with the micropipette along with the vesicle has diffused away 
enough not to affect the refracting index of the external solution significantly, and 
early enough that the amount of the permeable solute which has entered the vesicle can 
still be neglected. This procedure requires some experience and can not be automated. 
This does not present a problem with the experiment where membrane permeability 
due to peptide incorporation in the membrane is studied [10], as the permeability 
of membrane in the absence of peptides is very low and the experimentalist has 
plenty of time for the calibration. As the technique presented here does not rely 
on the properties of the halo profile, it provides a better alternative for measuring the 
permeability of membrane for small polar molecules. 

4. Conclusions 

We describe a technique for determining the permeability of a phospholipid membrane 
of micrometre-sized GUVs from a sequence of videomicrographs, which complements 
the existing methods of determining membrane permeability on either sub- micrometre- 
sized liposomes (SUV, LUV) or planar membranes. The technique relies on 
determining the radius of a spherical GUV with a great precision, and offers 
a substantial improvement over similar techniques which rely on quantifying the 
properties of the phase-contrast halo such as the height or width of its profile. When 
applied to a test set of data, the technique yielded consistent results. The obtained 
value is, however, « 40% lower than the results obtained on similar systems. Further 
investigations on a system comparable to the ones in the literature are needed to 
determine whether this discrepancy would still persist. 
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